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1 Introduction. 



The new HERA data JE] show the steep ^-dependence of the total cross section in the deep 
inelastic scattering ( DIS ) of virtual photon off a proton ( <J to t{.1*p) )• Approximately, 
c(7*p) oc x~ ' 2 at small x ( 10 -2 < x < ). Surprisingly, this energy rise holds at 

rather small photon virtualities ( Q 2 ~ 1 — 2GeV 2 ). At first sight it means that in HERA 
kinematic region we still have sufficiently diluted parton cascade and the parton-parton 
interaction which shall stop the increase of the parton density is still rather small. On 
the other hand, the probability of the parton - parton interaction 0] (« ) is equal to 



- = xG(x,Q)-^ = ^xG(x,Q 

where xG(x, Q 2 ) is the number of partons ( gluons) in the parton cascade, a is the cross 
section of parton-parton interaction and R 2 is the size of a hadron. The numerical factor in 
Eq. (fl|) has been evaluated by Mueller and Qiu || and has been confirmed in many further 
publications |4j]. Fig.|I] shows the contour plot for k using the GRV parameterization || for 
the gluon structure function and the value of R 2 = 5GeV~ 2 . We will argue a bit later that 
this value of R 2 follows directly from HERA measurement of the diffraction production 
of J/\l/ meson ||. One can see that k reaches k = 1 at HERA kinematic region, meaning 
shadowing corrections take place. Therefore, the situation looks very controversial. 

The goal of this letter is to derive the Froissart boundary for the deep inelastic structure 
function. This should clarify when the shadowing corrections to the deep inelastic process 
become important. Some attempts have been made to derive the geometrical limit ( the 
Froissart boundary) during the last three decades ( see Refs. || and lectures |§ for 
update review on the subject) assuming that a target is black for the dominant hadronic 
component in the wave function of the virtual photon. We derive the Froissart boundary 
for the deep inelastic processes assuming the GLAP evolution equation for the parton 
densities and the colour dipole picture of interaction proposed by A. Mueller |TT 



( see also Refs.[[L2|] [|13|] [|l4j] where many of aspects of the Mueller approach have been 
foreseen) . 



2 Unitarity constraint for DIS. 



2.1 s - channel unitarity ( general formulae). 

The unitarity constraint can easily be derived considering the DIS in the frame where a 
target is at rest. In this frame the virtual photon at high energy ( small x ) decays in 
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quark - antiquark ( qq ) pair long before the interaction with the target. The qq system 
traverses the target with fixed transverse distance r±_ between quark and antiquark [{12 



[|IIJ. Indeed, r± can vary by amount Ar± oc i?^, where Q denotes the energy of the 
qq pair or the virtual photon in the target rest frame, R is the size of the target, and the 
quark momentum is k± oc — ( see Fig.^|). Therefore 

k 

Ar ± oc R — < r± , (2) 
Qo 

which in terms of x has the following form: 

x « ■ (3 > 

The cross section for the DIS can be written in the form: 

<r{l*p) = £ dz J d 2 r ± q(z,r ± )a tot (zQo,rl)y*(z,r ± ) , (4) 

cross section for qq interaction with the target, z is the fraction of energy of the photon ( 
Qo ) carried by quark and r± is the transverse separation between quark and antiquark. 
\I/ is the wave function of qq - pair in the virtual photon. This wave function is well known 
[pTOfl |i~4| and \^t\ 2 for transverse polarized photon and for massless quarks is equal to 

G7YI A T ^ f 

I*t| 2 = ^^Z>{z> + (l-zf}Q>Kl(Qr ± ), 

where K\ is the modified Bessel function, Q 2 = z(l — z)Q 2 , Nf is the number of massless 
quarks and Zf is the fraction of the charge carried by the quark. 

The main contribution in Eq. ([|) (see Ref.[|T0[| for details ), which corresponds to the 
GLAP evolution, comes from the region Qr± < 1 and z{l — z) < qhpr < \- in this 
case the integral over z can be taken explicitly. Since z <C 1, it can be reduced to the 
integral 

/■i - - r°° ,— 8 

/ dzQ 2 K 2 (Qr ± ) -> 2Q 2 / z dz K 2 (</zQr±) = ——^ , (5) 

J J o LJ T i 



with Q r]_ > 4. Finally, we have 



At high energies ( low x ) we can restrict ourselves by summing only (as ln(l/x)) n contri- 
bution in each order of perturbative QCD ( so called leading log(l/x) approximation 
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(LL(x)A)). In the framework of the LL(x)A we can safely replace the argument of a in 
Eq. @ by x. Taking into account the relation between the cross section and F 2 structure 
function, namely cr to t(7*p) = 47T q? s ^(x, Q 2 ), the final formula has a form Jit] ]: 

The total cross section for qq scattering can be written as 

^.2 n. 

<7tot(x,—) = 2 J d 2 b ± Ima(x,r ± ,b ± ) , (8) 

where a is the amplitude for elastic scattering of qq in impact parameter (b±) space which 
is defined as 

a(x,r ± ,b ± ) = J- / <Pq ± e- i *^f(x 1 r ± ,t=-ql), (9) 
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where q± is the momentum transfer (see Fig.Q). In this representation 

a d = J d 2 b_L\a(x,r ± ,b_L)\ 2 . (10) 

The amplitude is normalized such that: 

^ = n\f(x,r ± ,t = -ql)\ 2 ; (11) 
a tot = 4<nlmf{x,r ± ,t = 0) . (12) 

The s - channel unitarity establishes the relationship between the elastic amplitude 
(a) and the contribution of all inelastic process ( Gi n (x,r^,b±) ) and has the form: 

2Ima(x,r_ L ,b±) = \a(x, r±, b±)\ 2 + G in (x,r±,b±) . (13) 

The solution of the unitarity constraint of Eq. ( [[3D is very simple if we assume that the 
elastic amplitude is predominantly imaginary at high energy. Indeed, one can check that 
the general solution of Eq. (|T3|) in this case has a form: 

a = i{l - e -|^x,6±)} ; (14) 

G m = {1 - e -n(*.'x,6x)} . (15) 

where fl(x, r±, b±) is the opacity function. Substituting Eq. ( |T4]) in Eq. (|j) and Eq. (0), 
we obtain 



3 



2.2 Properties of Q. 



The opacity f2 is an arbitrary real function which requires more detailed QCD calculations 
in order to be found ( see for example Refs. [0] [[| JT2| [|Kj]) and/or use the general property 
of analyticity and crossing symmetry ( see Refs. [|l8j |19| ). 

Let us recall what is known about Q: 

1. If Q <C 1 one can expand the exponent in Eq. (|14D and Eq. (|7]) can be reduced to 
a simple form: 

f^Q 2 ) = i^MjTf/ d% ± n( x , r l b± ) ■ (it) 

Differentiating over In Q 2 and comparing Eq. fll7|) with the GLAP evolution equations in 
the region of small x one can obtain for Q the following result ( see Refs. [0 JHJ \TE\ 
@ |7| for details): 

' d^fo ±b ± ) = 47T \ a ^ 2) x G(x, Q 2 ) , (18) 



/ 



where xG(x, Q ) is the gluon structure function of the proton. 

2. In the GLAP evolution equations the 6j_-dependence of the deep inelastic structure 
function can be factorized ( see Refs. 01 [OH in the form: 



47r 2 a s (Q^ 
3Q 2 



n = - ^ > xG(x,Q 2 )S(b 2 ± ) , (19) 
with the profile function S'(6 2 L ) equal 



S(bl) = JL J d 2 q± e-^F( q l) , (20) 

where F(t) is the two gluon form factor of the proton pictured in Fig.^a. Using, for 
example, the additive quark model (AQM) we can expect that this form factor is equal 
to the electromagnetic proton form factor ( see Fig.|3|b). Taking two different form of 
the proton form factor: the dipole ( Fdi P = ( 1 + R 2 q -2 j% p ) and exponential ( F exp = 

exp(—^R 2 qj_ ) ones, one can find two different profile functions, namely: 

S dip (b ± ) = ^2^^1(2^); (21) 
4 



and 



S exp (b ± ) = (22) 

7T iX z 



with normalization / d 2 b±S(b±) = 1. 

3. We can recover the eikonal (Glauber) model for the shadowing corrections (SC ) 
if we postulate Eq. (pT9|) for Q with the profile function S(b±) from Eq. (|2~ID or Eq. (p2]) 
for any values of b±. The physical meaning of this assumption is very simple: the final 
inelastic state is an uniform distribution that follows from the QCD evolution equations. 
In particular, we neglect the contribution of all diffraction dissociation processes to the 
inelastic final state that cannot be given as the decomposition of the qq wave function. 
For example, we neglect the so called "fan" diagrams (see Fig.|3]c) which give the dominant 
contributions at very large values of Q 2 and small x JJ. 

4. At large value of b± > bo±, Q falls down as Q oc e~ 2tlb± where fi is the mass of the 
lightest hadron ( pion ). Assuming that the DIS cross section cannot increase faster than 
(-) N , where power N = 1 comes from analyticity and crossing |18|] [19|, one can obtain 
the estimate for value of b± . Indeed, 



gives b ±0 = ^ ln± + 0(^j^] 



3 Froissart boundary for i^. 

Actually, Eq. (|23|) gives us the Froissart boundary for F2. Differentiating Eq. fllBf ) over 
InQ 2 we obtain (for A^ c = 3 and Nf = 3) 

dF 2 (x,Q 2 ) Q 2 r , 2 , -in (l ',y 



dlnQ 2 3tt 3 



J d 2 b ± {1 - e -h^M } (24) 



Using Eq. (|23[) we derive the estimate 



dF 2 (x,Q 2 ) Q 2 bj Q 2 N 2 2 1 2 1 

^Tng^ < = rW - Q x ' (25) 

where Q 2 is in GeV 2 . 

This boundary turns out to be well above all experimental observations. However, 
we can use a more detailed experimental information to obtain more restrictive estimate. 
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Indeed, the HERA data on diffractive photo and lepto production of vector mesons [50 
supports the idea that t-dependence of the DIS amplitude can be factorized out in the 
form of Eq. ( |T9"D or, in other words, as F 2 (x,Q 2 ;t) = F 2 (x,Q 2 ) ■ F(t) where the slope 
in t corresponds Eq. fl2T|) and/or Eq. (p2[). Using such form we can obtain directly from 
Eq. (J5JJ) the boundary 

dF 2 (x,Q 2 ;b ± ) _ Q 2 {1 _ e -|n(x,^), (2g) 



dlnQ 2 3tt 3 



Now, using 



F 2 (x,Q 2 ;b ± ) = F 2 (x,Q 2 )S(b ± ) < F 2 (x,Q 2 ) S(0) 
we derive from Eq. (p6|) 



dlnQ 2 3S(0)n 
which gives for the profile function from Eq. ( pT|) 



dF 2 (x,Q 2 ) < Q 2 ^ 



dF 2 (x,Q 2 ) Q 2 2 . 
«91nQ 2 < ^ > 
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while for Eq. (0) we have 



dF 2 (x,Q 2 ) 

dlnQ 2 < 3tt 2 ^ • 1 J) 

Taking R 2 = 10 GeV~ 2 which corresponds both to soft high energy phenomenology [21 
and the experimental data on diffractive lepto and photo production of vector mesons |20 
we are able to compare Eq. (f23) and Eq. (|29|) with the experimental data on 9F2 ^ X 



8F 2 (x,Q 2 ) 



dlnQ 2 



see Refs.p|] @ H). In Fig .| we plot the ratio R = ai *f , where FB is the Froissart 



boundary taking in the form of the Eq. (^Sj) or Eq. (^). For F 2 (x, Q 2 ) we used the GRV 
parametrization which fit the data quite well. One can see that the GRV parametrization 
reaches the unitary boundary ( R = 1 ) at Q 2 = Ql = 2 — 4 GeV 2 at HERA kinematic 
region. We can estimate the value of even more accurately using the parameterization 
of the experimental data on F 2 given in Ref . |[24|| , namely, 9F q^q2 — 0.364 log . Com- 

12 _ i nno^ 2 i„„ 0.074 



paring this parameterization with Eq. one obtains Q 2 = 1.092^ log™. Therefore, 
the value of Qq turns out to be pretty high at low x. This fact encourage us to search for 
a more microscopic approach for the parton - parton interaction in the parton cascade at 
moderate values of Q 2 ~ 2 GeV 2 . 
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4 Froissart boundary for the gluon structure func- 
tion. 



As has been pointed out by A. Mueller [[TO], the gluon structure function can be also 
written through the dipole GG - pair interaction with a target in a similar way as has 
been done for F%. Omitting all calculation that can be found in Refs. [|l(| 
derive 



25 



one can 



d 2 xG(x } Q 2 



<91nQ 2 «91n± 



7T" 



d 2 bl 



where the opacity Qgg f° r GG - dipole scattering has the same properties ( see section 2.2) 
as for qq- dipole scattering. The difference is that in the limit of small fi^G; ^gg — § ^qq 
for N c = 3. Repeating all arguments of section 3 one can obtain the Froissart boundary 
for xG(x, Q 2 ) in the form 



d 2 xG(x,Q 2 ) 
,9 In Q 2 ,9 In ± 



< 



2Q 2 



7T Z 



bio 



2Q 2 

4/i 2 7T 



in"! 

X 



2.5 Q 



2 In 2 - 

x 



(31) 



This boundary is much higher than the f-pTpg^ i i n a U current parameterizations of 
the gluon structure function |f26| || p7| . However, using the approach developed in the 
previous section one can obtain more restrictive estimates, namely 



d 2 xG(x,Q 2 ) Q 2 
,9 In Q 2 ,9 In ± < ^ 



I? 2 



and 



< 



2Q 2 



7T 



/? 2 



91nQ 2 aini 

and Eq. (^), respectively. In Fig.^] we plot '''''''' 



(32) 



(33) 



¥l for the 



for S(b±) from Eq. v _ i , rui-t J^[- \ — ;. m.-n ( m» <m .,<>. m J i^.-j m |mM- .,j J, 

GRV parameterization of the gluon structure function and compare them with Eq. (|3"^) 
and Eq. (^). One can see, that the gluon structure function reaches the unitarity limit 
( R — 1 ) at HERA kinematic region. 



In Ref.pl it has been shown that dxG \ x ^ ) 

11 11 am - 



3 from the experimental data on Fi- It 

means that the gluon structure function reaches the unitarity boundary at Q\ m ^ ~ 
3GeV 2 ( see Eq. 
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5 Parameter for the SC. 



To find out the parameter for the SC let us rewrite Eq. ( |30|) in the kinematic region where 
^gg < 1- One obtains 

d 2 xG(x,Q 2 ) 2Q 2 r 2 11 1 2 1 

d In In I = —J db±{ 2 n ° G ^ QS' b±) ~ 8 ^ g(X ' Q*' " (34) 

Substituting fi GG = xG(x, Q 2 ) 5(6 ± ) we have 

l^fg - tf) " W / *i } , (35) 

which can be rewritten in the form 

d 2 *G(x,Q 2 ) _ ^l xG (x O 2 ) • II - -) f36) 

din din ± ~ vr xG ^J 4 j ' W 

with 

-xG(x,Q 2 ) / rffei^ 2 ^) . (37) 



as^vr 3 



The above equation gives the same definition for k as Eq. (|l|) for exponential form 
of S(b±) ( see Eq. ( |22|) ). Using the new HERA data on photoproduction of J/\I/ meson 
|| we are able to estimate the value of i? 2 in the definition of k , recalling that R 2 is 
the size of the target only in the oversimplified eikonal ( Glauber ) model. To illustrate 
the point we picture in Fig.|| the process of J/\l/ photoproduction in the additive quark 
model (AQM ). We see that we have two processes with different slopes (B ) in t ( or 
in b\ ): the J/\l/ production without ( Fig.|a ) (B ei = 5GeV~ 2 ) and with ( Fig.|]b ) ( 
B in = 1.66 GeV~ 2 ) dissociation of the proton. The AQM gives us the simplest estimates 
for the resulting slope ( R 2 ) in Eq. (|l|) if we neglect any slope from the Pomeron - J/\l/ 
vertex in Fig.^, namely 

w = + *k } ~ ■ < 38 > 

This is a reason why we used R 2 = 5 GeV~ 2 in Fig.|I| to estimate the scale for the SC. 

In our estimates of the value of the deep inelastic structure functions at b± = ( see 
Eq.(27) ) we used an assumption that the SC does not change the value of R. To justify 
this assumption we plot in Fig. 7 the ^-dependence of the average b\ calculated in the 
Glauber (eikonal ) approach with R 2 = 5GeV~ 2 . One can see that < b\ > only weakly 
depends on x in the HERA kinematic region. 
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6 Summary. 



It has been presented the derivation of the Froissart boundary for the deep inelastic 
structure functions. 

The comparison of the Froissart boundary with HERA experimental data shows that 
both F 2 (x, Q 2 ) and xG(x, Q 2 ) hit the unitarity limit at Q 2 « 2 — AGeV 2 . This fact gives 
rise to a challenge for theoreticians to explain the matching between the deep inelastic 
scattering and real photoproduction process in the framework of QCD. 



We hope that this letter as well as Ref. [24| will stimulate the new round of the discus- 



sions on the theory of the shadowing corrections in the deep inelastic processes. We believe 
that the resolution of all difficulties could be found assuming that the SC has worked in 
the full in the gluon structure function and has been taken in the phenomenological initial 
gluon distribution in standard parameterizations || |26| [27j. However, much more work 
is needed to prove this. 

This work was supported in part by the U.S. Department of Energy, Division of High 
Energy Physics, Contract No. W-31-109-ENG-38 and by CNPq, CAPES and FINEP,Brazil. 
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Figure 2: DIS in the target rest frame. 
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Figure 3: Two gluon proton form factor in the additive quark model (a) and b)); c) The 
"fan" diagrams for DIS. 




QpGRV ( x q2\ 

Figure 4: The comparison of — 2 g\ n Q2 — given by Eq. ( p4j) with the Froissart boundary 



(Eq. (|28|) and Eq. (j29|) ) with the GRV parametrization. R 
Froissart boundary given by Eq. (|28| ) and Eq. (f29|). 



8F 2 (x,Q z ) 

FB 



where FB is the 
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R=d xG/dlnQ dln(1/x)/FB for McDonald profile. 

R 2 = 10 Gev 2 . 



Q = 1 . Gev 

= 2.5 Gev 2 

=5. Gev 2 

= 1 0. Gev 2 

= 20. Gev 2 




R=d 2 xG /dlnQ 2 dln(1/x)/FB for exponential profile. 



Q = 1. Gev 
= 2.5 Gev 2 
• =5. Gev 2 
■ =10. Gev 2 
= 20. Gev 2 




y=ln(1/x) 



10.0 15.0 20.0 

y=ln(1/x) 



Figure 5: The comparison of 



given by Eq. (gD and Eq. fl33|). R 
by Eq. (||) and Eq. (H) 



9 Q^Q2gi n ( f ' given by Eq. fl3"0|) with the Froissart boundary 

d 2 xG(x,Q 2 )^ 

a in y in - ^ w j iere j g p ro j ssar t boundary given 





Figure 6: The J/^f production without a) and with b) dissociation of the proton. 
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Figure 7: The x - dependence of < b]_ > in the Glauber (eikonal) approach with R 2 = 
5GeV~ 2 . 
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